Peptide size dependent active transport in the proteasome 
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We investigate the transport of proteins inside the proteasome and propose an active transport 
mechanism based on a spatially asymmetric interaction potential of peptide chains. The transport 
is driven by fluctuations which are always present in such systems. We compute the peptide-size 
dependent transport rate which is essential for the functioning of the proteasome. In agreement 
with recent experiments, varying temperature changes the transport mechanism qualitatively. 

PACS numbers: 05.40.-a,05.40.Ca,05.60.-k 



Eukaryotic mammalian proteasomes are fascinating 
molecular machines, intended to digest tagged with ubiq- 
uitin proteins [lj . The proteins are cleaved into peptides 
which can be further used for the synthesis of new pro- 
teins. The proteasome's function has been directly linked 
to a wide range of diseases, such as cancer, neurologi- 
cal diseases, and diseases of the immune defense system. 
Starting from the discovery of their activity in the late 
1970's, presently proteasomes are in the focus of cur- 
rent molecular biology, see and refs. therein. By now 
there are several models to explain theproteasome's ac- 
tivity H, cleavage mechanisms [J EJ and to predict 
the cleavage results 0, however, basic principles of the 
proteasome operation mechanisms are still poorly under- 
stood, mainly due to the lack of experimental results. 

In this Letter, we focus on understanding of protein 
translocation inside the proteasome, leaving the mecha- 
nisms of cleavage, targeting, etc. beyond the scope. The 
main question is: given proteasomes are highly complex 
pipe-like structures of tenthousands atoms of almost per- 
fect left-right symmetry with respect to the axis of the 
pipe 0, thus, operating bi-directional. To be cleaved, 
a protein has to enter the proteasome at one side, pass 
the active sites (where the cleavage occurs) , which makes 
about 1/3 of the total length of the pipe. Then the 
cleaved peptides have to pass all the way through the pipe 
to finally reappear at the other side of the proteasome. 
Although there is a number of examples where protein 
transport in cells occurs due to diffusion, i.e. Brownian 
motion Q, diffusive transport may be excluded as the 
main mechanism for translocation in the proteasome be- 
cause of the enormous cargo [ljj. Also other proposed 
transport mechanisms, such as the power stroke model of 
protein translocation do not suffice to explain transloca- 
tion Therefore, the question arises how the protein's 
motion inside the proteasome is driven? 

Since proteasomes are large multi-subunit structures 
consisting of proteins, the mechanism of the protein 
transport is directly related to protein-protein interac- 
tion. In 12] it has been noticed that if attachment 
and detachment rates are specified asymmetrically, the 
protein-protein binding interaction acts as a ratchet. 
Following this argumentation, active protein transport, 
based on the mechanism of a molecular ratchet, has been 



studied for transport through membranes |13| and has 
been also discussed as a mechanism for cytosolic destruc- 
tion by the proteasome |14| . Moreover, maximum likeli- 
hood tests have shown that other models, e.g. the power 
stroke model of protein translocation, do not lead to bet- 
ter agreement with the experiment than the Brownian 
ratchet model Noteworthy, in these ratchet effects 
transport is possible only in a certain temperature inter- 
val, and stochasticity, intrinsically present due to fluctu- 
ations in any biochemical reaction |l5l llf^ | , provides the 
driving mechanism. 

Here we propose a model for active protein transloca- 
tion in the proteasome to explain the peptide size depen- 
dence of the transport velocity as well as its tempera- 
ture dependence which possibly explains the mechanism 
of temperature reaction or heat shock response, regulat- 
ing the proteasome activity in the case of some diseases 
|17llialial2fl| . The results describe a system size ratchet 
effect that is related to similar effects which have been de- 
scribed recently for other noise-induced phenomena [2l| . 
It has been shown in vitro (e.g. 0) that already the 20S 
proteasome reveals full functionality, thus, for studying 
the transport we disregard the enzymatic action of the 
proteasome caps. 

Model. Following, e.g., we assume that the protein- 
proteasome interaction is characterized by a spatially pe- 
riodic asymmetric potential, motivated by asymmetric 
C-N binding between amino acids that enter the protea- 
some always in N — > C direction, i.e., with N head on 
@ (see Fig. P). The folding structure of the protein is 
not relevant here since due to the action of the regula- 
tory complexes the protein enters the proteasome unfold 
|22j . Moreover, even fold proteins tends to form periodic 
structures to maximize their amphiphilicity p3| . 

As an abstraction, here we assume that the protein- 
proteasome interaction potential U{x) is periodic with 
the period L. In reality there is a basic asymmetry, 
namely the C — N asymmetry of the protein (or pep- 
tide) backbone, that is superposed by a non-periodic (in 
our sense irregular) part that is attributed to the amino- 
acid-specific residues [24] . The basic structure of a pep- 
tide and its model as an asymmetric ratchet potential is 
sketched in Fig. ^ The commensurability of the pro- 
tein and the proteasome is supported by the fact that 
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FIG. 1: The regular part of the protein-proteasome interac- 
tion potential is modeled by a spatially asymmetric periodic 
potential, based on the C — N asymmetry of the peptide back- 
bone. The Ri mark the amino-acid-specific residues. The 
lower part shows the equipotential surfaces of a simple pep- 
tide which is clearly asymmetric (see also |24j|1. The H and 
O atoms are marked for orientation. 



both biological macro-molecules consist of the same basic 
structures, namely chains of amino-acid. That commen- 
surability has been claimed before, and was the basis of 
automaton-like models of proteasome digestion processes 
[f| . The detailed form of the asymmetric periodic inter- 
action potential is of less importance for this qualitative 
study. Here, we assume a saw-tooth potential as drawn 
in Fig. The angles are smoothed, i.e., dU/dx can be 
computed in each point (for details see |25f1. 

The proteasome acts upon the protein by a certain 
number of equidistant interaction centers. The dynam- 
ics of the protein inside the proteasome is, hence, gov- 
erned by N interactions centers, where N is the num- 
ber of protein elements (amino acids or multiplicatives 
of it). There appear the following forces: potential force 
(protein-proteasome interaction) —NdU(x)/dx, fluctua- 
tions with collective NF(t) and individual components 
fi(t) + ... + /jv(t), and protein friction forces Nf3x [l2| . 
where x is the coordinate of the protein with respect to 
the proteasome and f3 is the coefficient of friction. Due to 
the small mass of all protein particles, moving in the liq- 
uid cytosol, the motion occurs in the overdamped realm 
|26j. hence we neglect inertia forces. Note that transport 
is possible only in the case of nonequilibrium fluctua- 
tions. In the simplified case, when fluctuations can be 
represented by a sum of a collective periodic force and 
individual for every protein residue thermal noise, the 
model is analytically tractable, predicting the velocity 
dependence on the peptide size. Numerically, we inves- 
tigate also different kinds of fluctuations in order to be 
closer to reality. Normalizing all forces by friction and 
taking = 1, the translocation of a protein in the pro- 
teasome is governed then by 
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+ F(t) + -(Mt) + ... + f N (t)). (1) 



FIG. 2: From left to right; Top: protein moves into the pro- 
teasome while increasing the number TV of interaction centers. 
Bottom: after cleavage N abruptly decreases for the resulting 
peptides. When leaving the proteasome, N decreases as well. 



When the protein chain enters the proteasome, the num- 
ber of interaction centers N is increased (Fig. EJ). After 
cleavage, N is abruptly decreased for the cleavage prod- 
ucts. We believe that this crucially changes the trans- 
port velocity, providing the cleavage product leaving the 
proteasome. Eq.0 models the motion for both the ini- 
tial protein and the peptides after cleavage. We consider 
three cases of fluctuations, regulating their motion: 

Case 1: We assume collective oscillations of the pep- 
tide elements jHH^, e.g., F(t) = Acos(ajt), where A 
and u! stand for the amplitude and frequency of these 
oscillations. Additionally, each interaction center under- 
goes local thermal fluctuations, represented by mutually 
uncorrelated white noise of intensity a 2 : fi(t) = £j(t), 
where (&(t)&(f)) = a 2 5(t - t')8 l] . In this case the 
stochastic term in Eq. Q is white noise of intensity 
a 2 /N. The Fokker-Planck equation for the peptide co- 
ordinate probability distribution w(x,t) associated with 
Eq. © is 
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which may be solved in quasi-stationary adiabatic ap- 
proximation dw/dt — I23. We obtain 
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w(x,F)=-G(F), (2) 



where G(F) is the probability flux. For any periodic 
potential U (x) the quasi-stationary solution of Eq. is 
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where C(F(t)) and G(F(t)) are unknown functions of t. 
Using the periodicity condition w(0,t) = w(L,t) and the 
normalization of w(x, t) we get G(F). If the amplitude A 
meets the condition LA <C a 2 /N, one can expand G(F) 
and obtain 
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FIG. 3: left: Mean velocity versus noise intensity for N = 
1, 2, ..7. right: Mean velocity vs N for different temperatures: 
a 2 = 0.6(curve 1), 2, and 3 (curves 2 and 3). The inset plot 
shows the depth of the protein penetration for velocity 0.0005. 



with the expansion coefficients Gni = L/ (h^ho), 
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Substituting Eq. © into (x) — L G(x,t)dx, where (•) 
denotes time averaging, we obtain the average protein 
transport velocity, as a function of the noise intensity a 2 
and the peptide size N 
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The dependence of the transport velocity on the noise 
intensity, i.e. temperature, is shown in Fig. [3] (left) for 
different system size N . This plot shows that protein 
transport is possible only in a certain temperature in- 
terval. Absence of noise leads to vanishing transport as 
well as too high temperature also decreases the veloc- 
ity. For given temperature, the peptide size N strongly 
influences the transport velocity, i.e., with increasing N 
the dependence is shifted to the right and becomes more 
wide (curves 1-3). This dependence differs qualitatively 
for small and large temperature. For rather small tem- 
perature the velocity decays with increasing N (curve 1 
in Fig. [3] (right)). Consequently, when a protein enters 
the proteasome, it will initially move with rather large 
velocity. Then, due to the increase of N the velocity 
decreases, and the motion of the protein virtually termi- 
nates. On this stage the probability of cleavage is much 
higher than that of translocation, thus, with large proba- 
bility the protein will be cleaved there. After cleavage, N 
is significantly decreased (Fig. [21 (right)), and following 
the velocity dependence, both fragments move rapidly to 
leave the proteasome. Noteworthy, in agreement with 
our model, the blocking of uncleaved proteins inside the 
proteasome has been observed also experimentally |30| . 



Surprisingly, for large temperature the velocity de- 
pends qualitatively different on the peptide size (curves 
2,3 in Fig. |3 (right)), i.e., there is an optimal peptide size 
that corresponds to maximal velocity. In this case the 
protein moves accelerated while entering the proteasome 
until the optimal velocity is achieved. Then the velocity 
decreases, and again it virtually stucks. After cleavage 
the behavior is quite different: if the cleaved peptide is 
close to the optimal size, it leaves the proteasome rapidly. 
If the cleaved fragment is too small or too large, it moves 
with small velocity, thus blocks the proteasome leading 
to low efficiency or malfunction of the proteasome. Note 
that for the protein that enters the proteasome the initial 
small N increases, and the velocity will increase, whereas 
for the cleaved fragment of smaller N the velocity re- 
mains constant. There is experimental evidence |30| that 
certain proteins cannot be cleaved by the proteasome, 
i.e., they are inert with respect to its active sites. For 
such proteins we can determine the penetration depth, 
defined as the depth N at which the propagation veloc- 
ity falls below a predefined very small threshold. The 
inset of Fig. (right) shows that the penetration depth 
increases linearly with temperature. So far, there are no 
experimental data to verify this prediction. 

Consequently, small temperature leads to preferably 
small size of the cleaved fragments, whereas larger tem- 
perature leads to domination of some specified peptide 
size in the proteasome output. As a hypothesis, we be- 
lieve that the qualitatively different proteasome behavior 
with respect temperature may be responsible for temper- 
ature reaction or heat shock response 0, 0, |2(j , pro- 
viding immune defense in the case of some diseases. For 
large temperature (Fig. [3] (right) curve 2,3) the protein is 
accelerated while entering the proteasome, hence, leading 
to more efficient operation of the proteasome. 

To confirm these results we have performed computer 
simulations [3l| whose results agree with our analytical 
findings (Fig. |3]|. Some quantitative deviations result 
from taking into account only a limited number of terms 
in the expansion Eq.JSJ). 

Case 2: Since collective fluctuations of the protein or 
peptide elements, represented by a periodic in time force, 
can be used only under certain assumption, next we con- 
sider the system in the absence of deterministic periodic 
forces, but with colored noise F(t) — ((t), where ((t) is 
harmonic noise with some dominating frequency (due to 
the molecular spectrum and oscillations of the peptide 
chain H^H^)- We obtain this noise as the solution of 



C + 2<5C + ^ 2 C = £', 



(5) 



where £' is Gaussian noise {£'(t)£'(t')) — cr 2 8{t — t'), and 
parameters S = 0.01, uj = 0.1, a 2 = 10~ 4 regulating the 
width, the dominating frequency of the power spectrum, 
and the intensity of ((t). The results show qualitatively 
the same behavior as in the case of the periodic force 
(Fig. 0] (b,e)). Hence, if exclusively noisy fluctuations 
are present in the system, the generic behavior described 
above, persists and leads to the same temperature and 
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FIG. 4: a,b,c: Mean velocity versus noise intensity for differ- 
ent N =1,2, ..5. Here A = 1.15, w = 0.1. d,e,f: Mean velocity 
vs. N for different temperatures: d) a 2 = 0.3 (1), 2 (2); e) 
a 2 = 0.3 (1), 1 (2), and 2 (3); f)cr 2 = 0.5 (1) and 1.5 (2). 
a,d: with periodic force, b,e: with colored noise, c,f: without 
global fluctuations. 



peptide size dependences. 

Case 3: Next we consider the case of no collective fluc- 
tuations, F(t) = 0, i.e., there is no eigen collective mode 
in the peptide oscillations. The local noise of each inter- 
action center is presented now by the sum of uncorrelated 
colored and white noise components: fi(t) = Q(t)+£i(t). 
The presence of the colored noise is necessary to generate 
nonequilibrium fluctuations, that is a necessary condition 
for directed transport. In this case the motion of the pro- 
tein or peptide is described by 



dU(x) 
dx 



1 

N 



(Ci 



+ Cn + fl 



•60 
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with Q(t) and £j(t) as above. In general, even in this 
case the behavior is similar, but increase of the peptide 
size leads not only to a shift of the velocity dependence 
towards large noise intensity, but also to a decrease of its 
maximal values (Fig. 01(c)). This is due to the fact that 



nonequilibrium fluctuations originating from different in- 
teraction centers act occasionally in converse direction. 
For small temperature (Fig. 01 (f), curve 1) the trans- 
port decreases with increasing peptide size N, whereas 
for large temperature the rate depends nonmonotonously 
on N (curve 2). The optimal velocity in the latter case 
is smaller, as for the case of collective fluctuations. 

In summary, we investigated a model for active protein 
transport inside the proteasome analytically and numer- 
ically. The size of the peptide as well as temperature in- 
fluence the transport rate significantly in agreement with 
experimental results. Following our predictions, protein 
transport occurs only in the certain interval of temper- 
ature. The size dependence of the transport velocity 
leads to a preferred fragment (peptid) size since for such 
sizes the probability of cleavage exceeds the probability 
of further translocation significantly. Under certain con- 
ditions, uncleaved proteins may get stuck inside the pro- 
teasome, thus leading to decrease of its efficiency which is 
in good agreement with experimental observations |30| . 
Calculating the average velocity, we have considered only 
the periodic constituent of the protein-proteasome poten- 
tial. Slight inhomogeneity (nonperiodicity) of this po- 
tential will lead to fluctuations of the velocity around its 
average leaving the qualitative behavior of dependencies 
on the peptide size intact. 

Varying temperature, the model predicts qualitatively 
different transport regimes, thus, possibly explaining the 
mechanism of temperature reaction or heat shock re- 
sponse as it is observed for certain diseases. For larger 
temperature, the protein is accelerated while entering the 
proteasome, moreover, there is a preferred fragment size, 
that should dominate the proteasome output. We em- 
phasize that our model is based only on two assumptions: 
spatial asymmetry of the protein chain and presence of 
nonequilibrium fluctuations, which seem to be certainly 
fulfilled in real proteasome machines. 

We thank C. Frommel, H. Holzhutter, and R. Preifiner 
for helpful discussion. 
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